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Abstract 

The renormahzation group method enables one to improve the properties of the QCD per- 
turbative power series in the uhraviolet region. However, it ultimately leads to the unphysical 
singularities of observables in the infrared domain. The Analytic Perturbation Theory con- 
stitutes the next step of the improvement of perturbative expansions. Specifically, it involves 



<^ \ additional analyticity requirement which is based on the causality principle and implemented in 

' the Kallen-Lehmann and Jost-Lehmann representations. Eventually, this approach eliminates 

^ ■ spurious singularities of the perturbative power series and enhances the stability of the latter 

, with respect to both higher loop corrections and the choice of the renormahzation scheme. 

The paper contains an overview of the basic stages of the development of the Analytic Per- 
turbation Theory in QCD, including its recent applications to the description of hadronic processes. 
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Preamble 



The method of the Analytic Perturbation Theory (APT) resolves the problem of un- 
physical (or ghost) singularities of both the invariant charge of Quantum Chromodynamics 
(QCD) and the matrix elements of the strong interaction processes. This difficulty (known 
also as the problem of Moscow zero or Landau pole) first appeared in Quantum Electrody- 
namics (QED) in the mid-50s of the last century. It played a certain dramatic role in the 
development of Quantum Field Theory (QFT). 

In the late-50s Bogoliubov, Logunov, and Shirkov suggested [H resolving this problem by 
merging the renormalization group (RG) method with the Kallen-Lehmann representation, 
which implies the analyticity in the complex Q^-variable. The method of APT in QCD is 
based on the ideas of Ref . p| . 

The development of the APT over the last decade has revealed a number of new prin- 
cipal features of the analytic approach. Specifically, in addition to the resolution of the 
problem of unphysical singularities, the APT leads to the nonpower functional expansion 
for QCD observables. The latter possesses an astonishing (in comparison with the pertur- 
bative power series) stability with respect to both higher loop corrections and the choice of 
the renormalization prescription. 

I. INTRODUCTION 

The renormalization group method, which was devised in the mid-50s in Ref. 0] (see also 
paper [3] and the respective chapter of book [3]), is an inherent part of the contemporary 
QFT calculations. This method becomes especially useful for singular solutions, when the 
type of the singularity is affected by perturbative expansion. Besides, the RG method is 
important for strong interactions, e.g., for QCD. 

The QCD description of the majority of hadronic processes requires the use of the RG 
method. At the same time, the straightforward solutions of the RG equations suffer from 
the spurious singularities. The one-loop QCD invariant charge a^iQ'^) possesses the ghost 
pole at Q"^ = A^, see Eq. fl2.ip below. Higher loop corrections just give rise to additional 
singularities of the cut type and do not eliminate this problem. The existence of such 
singularities contradicts the general principles of the local QFT 0], 

A solution to the problem of the unphysical singularities of the invariant charge was pro- 
posed in Ref. Specifically, this can be achieved by merging the RG method with definite 
properties of the analyticity in Q^-variable. The latter follows from the Kallen-Lehmann 
spectral representation for the transverse Lorentz-invariant amplitude of the dressed photon 
or boson propagator 



that reflects the basic principles of the local QFT. 

In QED, the square of the electron effective charge q;((5^), which was first introduced by 
Dirac [6], is proportional to the transverse amplitude of the dressed photon propagator. The 
latter satisfies the spectral Kallen representation (II. ip . which implies the analyticity in the 
complex Q^-plane with the cut along the negative semiaxis of real Q^. The function a{Q'^) is 




(1.1) 
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also called the invariant charge or running coupling constant^. In accordance with paper 
the analytic invariant charge can be reconstructed by making use of the representation (11.11) . 
the spectral density p{a) being defined as the discontinuity of the perturbative QCD invariant 
charge across the physical cut along the real negative semiaxis Re < 0. 

The QED analytic invariant charge d, elaborated in Ref. [l|, possesses the following 
important properties: 

- it has no ghost pole; 

- as a function of a it has an essential singularity of the form exp(— Stt/q;) at the origin; 

- for real and positive a it admits a power expansion, that coincides with the perturbative 
one; 

- a has finite ultraviolet limiting value Svr, which is independent of the experimental 
value a ~ 1/137. 

In the mid-90s this idea was employed in QCD in Refs. 0], Afterwards, the method 
developed therein was named APT. In QCD, the synthesis of the renormalization invariance 
with the Q^-analyticity has revealed a number of important features of the analytic invariant 
charge 0], [8|. In particular, a^iQ"^) has the wnroersa/ infrared (IR) stable point. Its value 
«e(Q^ = 0) is determined by the one-loop /5-function coefficient Pq. The value a-E{Q^ = 0) 
is a scheme- independent quantity, since it is not altered by multi-loop corrections. The IR 
limiting value does not depend on the scale parameter A, which can be evaluated by making 
use of experimental data. The set of curves aE^Q^), corresponding to different values of 
A, forms a bundle with the common point aE(0) = l//?o- Therefore, the imposition of the 
analyticity requirement essentially modifies the IR behavior of the analytic invariant charge. 

Another important feature of the analytic approach is that it enables one to define the 
invariant coupling q;m(s) in the timelike (Minkowskian) domain in a self-consistent way j^. 
Usually, contemporary QFT calculations involve explicit expressions for observables and 
other auxiliary RG-invariant (or covariant) quantities, being expressed in terms of the invari- 
ant charge. To achieve it, the quantity at hand has first to be represented in an appropriate 
form. For example, only the quantity defined in the Euclidean region can be expressed in 
terms of as{Q^)- 

Meanwhile, a number of observables (e.g., the effective cross-sections and quantities re- 
lated to inclusive decays) are functions of the timelike argument s = —Q^, with s being 
the center-of-mass energy squared. However, in the framework of the RG method, one 
cannot straightforwardly substitute the spacelike (Euclidean) argument with the timelike 
(Minkowskian) one. Indeed, in accordance with (11.11) the amplitudes of propagators, sim- 
ilarly to the transverse photon amplitude d{Q'^), and the relevant matrix elements acquire 
complex values for real negative (i.e., for the timelike argument). 

Interrelations between the RG-invariant quantities in the Euclidean and Minkowskian 
domains can only be established by making use of the linear integral transformations. The 
analytic properties of the invariant charge, which are violated by perturbation theory (PT), 
and can be recovered within the analytic approach afterwards, play a crucial role here. 

Both, the Minkowskian ayiis) and Euclidean a-^iQ"^) charges possess the common IR 
stable point a^^s = 0) = a-^iQ"^ = 0) = l//9o- In the ultraviolet (UV) region these couplings 
also have the same asymptotic behavior. However, functions a-^^Q"^) and ayi^s) cannot be 
identical in the entire energy range due to certain general arguments [lo|. 





We shall not use the term "running coupling constant" here due to its semantic nonsense. 
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A key feature of APT is the transformation of series in powers of as for observables into 
the nonpower functional expansions. The latter displays both a milder dependence on the 
choice of the renormalization scheme and an improved numerical convergence. 

Thus, one arrives at the self-consistent method of description of observables. This ap- 
proach possesses the renormalization invariance and is free of unphysical singularities and 
related difficulties. 



II. ANALYTIC PERTURBATION THEORY 



A. Singularities of the effective charge 



At the one-loop level, the RG summation of the ultraviolet logarithms leads to the singular 
expression for the QCD running coupling ( "invariant charge" ) 



a. 
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where Uf denotes the number of active flavors. The scale parameter is deflned here in the 

1 

well-known way, namely A = /ie ^"t^i^o . 

The function (12. ip has the IR unphysical singularity at = A^. This problem cannot 
be solved HI] by taking into account higher- loop corrections. Indeed, the latter just modify 
the type of singularities, but do not eliminate them. For example, the common form of the 
two-loop perturbative coupling 
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(2.2) 



possesses an additional unphysical cut due to the double- log dependence on Q^. 

A similar situation takes place in QED as well. It is worthwhile to note here that in this 
latter case unphysical singularities correspond to huge energy scales which have no physical 
meaning. At the same time, in the case of QCD the value of the parameter A is about 
several hundred MeV, that is within the physically-accessible range of energies. 



B. Analytic approach 

In general, the renormalized perturbative expansion can be further modifled in a certain 
way. For example, the Adler function D{Q'^) is representable by the double series in powers 
of the running coupling at a normalization scale /i^ and in powers of ln((5^//i^): 

^p*(f '"^)=E<i:^-^i-'^(f)- (2-3) 

In the UV domain this expression is ill-deflned due to a large value of logs^. 



^ This series is also meaningless in the IR domain due to the same reason. Besides, one might expect the 
factorial growth of the expansion coefficients dn,k at large n. 
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The RG method improves the properties of the perturbative expansion in the UV re- 
gion by accumulating "large logs" into the invariant charge OisiQ"^) = q;s(Q^//^^, In 
particular, the RG-improved Adler D-function becomes a function of dg only, and can be 
represented by the power series 

DMQ') = D{as{Q')) = Y,<oasiQ'y, PT = pt + RG. (2.4) 

n 

Unlike Eq. (12. 3p . the terms of Eq. (12. 4p are //-independent, and the expansion parameter 
ctsiQ"^) vanishes when — >■ oo, in agreement with the asymptotic freedom. 

At the same time, expansion (12.41) formally remains ill-defined in the IR domain due to 
unphysical singularities of the expansion parameter as{Q^)- Thus, on the one hand, the 
RG improvement of the perturbative series (12. 3p results in a crucial physical property of the 
asymptotic freedom. On the other hand, an important property of D-function, namely, the 
analyticity in the complex Q^-plane with the cut along the negative semiaxis of real Q^, is 
lost in Eq. (1^ . 

The APT method constitutes the next step in the improvement of the perturbative ex- 
pansion. This method employs the principle of renormalization invariance together with a 
fundamental principle of causality which is realized in the form of the Kallen-Lehmann inte- 
gral representation (II. ip . In the framework of the APT, the D-function can be represented 
by the nonpower functional expansion 

^pt(Q') = dna:{Q') ^ /^apt(Q') = dnAniQ'). (2.5) 

n n 

The Euclidean functions An{Q^) satisfy the Kallen-Lehmann representation 

A(g^) = - r^^^T?^, (2.6) 

TT Jo (T + Q^ 

with the spectral function being defined as the discontinuity of the respective power of the 
invariant charge across the physical cut: pfc(o") = lin.a^{—a — it). The first-order function 
AiiQ"^) corresponds to the analytic running coupling 

«e(Q') = - r da p{a) ^ p,ia). (2.7) 

The representation (12.60 determines the analytic properties of the functions AniQ"^) in the 
complex Q^-plane. Specifically, these functions (and, therefore, the series (12.50 ). are analytic 
functions in the complex Q^-plane with the cut along the negative semiaxis of real Q^. The 
investigation of the properties of these functions has revealed that the resolution of the 
ghost pole problem by making use of the APT method eventually leads to the IR stability 
with respect to higher loop corrections. The results of this investigation are summarized 
in Table [T] which elucidates the stages of the evolution of the perturbative results: PT — »• 
PT+RG ^ APT. 

The last row of the Table contains two sets of functions, namely, An{Q'^) and 2l„(s). 
The latter appear in the description of the processes depending on the timelike (i.e., 
Minkowskian) momenta. These functions naturally emerge in the study of the Drell func- 
tion R{s) which is the ratio of the inclusive hadronic cross-section of the process of the 
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TABLE I: Properties of various approximations 



Method 


Type of approximation 


Properties 






UV 


IR 


Analyticity 


PT 


Double set in powers of 










and InQ^/^^ 






+ 


PT + RG 


Power series in 










invariant charge q;s(Q^) 


+ 






APT = PT + RG + 


Nonpower functional expansions 








+ analyticity 


in AkiQ'^) and 2lfc(s) 


+ 


+ 


+ 



e"'"e~-annihilation to the leptonic one. Here the timelike argument s is the center-of-mass 
energy squared. The perturbative approximation of R{s) by the power series in the invariant 
charge, Rpt{s) = X]„?"na"(s), violates the relation between the functions Di^Q"^) and R{s) 

In the framework of the APT, the function R{s) takes the form of the functional expansion 

i?APT = 5^C?n2l„(s), (2.9) 

n 

with the coefficients dn being identical to those of Eq. fl2.5p . The relation between Ak{Q'^) 
and 

is similar to the relation (12.81) between D{Q'^) and R{s). The functions 21^ can also be 
expressed^ in terms of the spectral function Pfc(cr) [9|]: 

Ms) = - / —Pk{a). (2.11) 
In the leading order the Euclidean running coupling reads Q 



A2 

+ 



/5oLln(gVA2) A2-g2j- 



(2.12) 



Figure [H depicts its behavior for A = 200 MeV and A = 400 MeV. For comparison, the 
corresponding perturbative curves are also plotted therein. The enhanced stability of the 



Schwinger argued that in QED the RG /3-function is proportional to the spectral function of the 
photon propagator (i.e., invariant charge). However, this hypothesis is violated beyond the two-loop level. 
In the framework of the APT, the Schwinger's assumption turns out to be realized in a "hybrid form" for 
the /9- function corresponding to the Minkowskian charge (|2.13p . Indeed, the logarithmic derivative of the 
latter is proportional to the spectral function of the Euclidean charge, see Eq. (|2.1ip . 
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FIG. 1: Euclidean and perturbative Og charges. 
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FIG. 2: Higher loop stability of the Euclidean charge oe- 



APT expressions (in comparison with the perturbative ones) is demonstrated in Fig. [2], where 
the one-loop and two-loop functions are shown. The three-loop analytic running coupling 
practically coincides with the two- loop one (within the accuracy of 1-2%). 

Thus, contrary to the RG-improved perturbation theory, the analyticity, which emerges 
from the causality, leads to the stabilization of the behavior of the invariant charge in the IR 
domain. A key property of the approach at hand is that all the expansion functions assume 
the universal value at = which eventually results in the stabilization mentioned above. 
At the same time, the stability in the UV domain (starting from the two-loop level) is due 
to the asymptotic freedom. 

In the framework of the APT the invariant charge in the Minkowskian region can be 




defined 13|] in a self-consistent way as tlie first of tlie functions (12.110 : 



71 J ^ a 

In tlie leading order this function reads 



Q^M ( -5 ) = ~ — arccos 



1 d(7 

amis) = - — p(cr), p(o-) = pi(o-). (2.13) 



1 TT , S , , 

arctan— , L = ln— . (2.14) 



L>o /^ovr L' 

The Euclidean ( 12. 7p and Minkowskian (I2.13p charges share the same IR limiting value 

aE(0) = aM(0) = (2.15) 
Po 

which is independent of the loop level. 



C. Higher APT expansion functions 

These functions are necessary for the analysis of observables by making use of Eqs. (12. 9p 
and (12.51) . They satisfy the recurrent relations 

n>l ^ n>l 

which can be employed for their iterative definitions. To achieve it, one has to explicitly 
solve these relations by making use of additional assumptions of the form Pe>em = and 
Ak>K+i = 0. At the one-loop level (/?£>i = 0) the APT formulae have a simple and elegant 
form. In this case, proceeding from the first functions (12.121) and (I2.14p . one can show by 
making use of relations (I2.16P that 



2(e'-l)3y'' ^3 (^)- ^3 (^2 + ^2)2' ^-^\a2 



(2.17) 



The two-loop level is more complicated technically. The point is that the exact solution 
for as is expressed in terms of the special Lambert function here, which leads to cumbersome 
explicit expressions for Ak and [14]. 

Nonetheless, all APT functions obey the following important properties at any loop level: 

• Unphysical singularities are absent, no additional parameters being introduced. 

• Higher functions, (I2.17P etc., are not equal to powers of the first ones (I2.12p . (12.141) . 
They oscillate in the vicinity of ~ A^ and vanish in the IR limit, see Fig. ^ below. At 
the same time, these functions tend to the powers of ag in the UV asymptotic. 

• The expansions of observables in powers of as((5^) (for the Euclidean case) and in 
powers of as{s) (for the Minkowskian case) are replaced by the expansions over the sets of 
{Ak{Q^)} and {2lfc(s)}, respectively. The latter expansions display a faster convergence 
with respect to that of the perturbative case. 
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D. Analyticity in the Og-plane 



As it has been noted above, the expression for the Euchdean couphng a^\Q'^ / ,a^j) 
contains the nonanalytic term of the form exp[— 1/ ag] . The latter corresponds to the essential 
singularity at the origin of the complex ctg-plane. 

More than half a century ago, proceeding from a general reasoning, Dyson argued [l^ that 
such singularity appears in QED inevitably. The explicit form of this singularity coincides 
with both the term determined in Ref. ^] by merging the renormalization invariance with 
the causality condition and with the results of Ref. [13] obtained by making use of the 
functional saddle-point method. 

It is worth noting also that the logical inevitability of the nonpower type of the functional 
APT-expansions was discussed in detail in Refs. [ISJ . 

The conversion of the common QCD running coupling as{L) into the Euclidean a-E{L) 
or Minkowskian ayiiL) one is equivalent to the introduction of a new expansion parameter. 
It is worth considering two examples 

a^wf(a) = 4rarccos^ ^ a ^ = « + ^(1 - e^/^^o^V^ (2.18) 

The first one is similar to the choice of another renormalization scheme, since the function 
w^{a) can be expanded in powers of a. The other one transforms into identity in the weak 
coupling limit, since its second term e"^/*^^""^ does not contribute to the expansion in powers 
of a at a — >■ 0. 

At the same time, the conversions fl2.18p give rise to the transformations of the invari- 
ant charge a"s\L) and ai^\L) A^i\l) which are equivalent to Eqs. fl2.12p 
and flCTj) . 

The function wf (ctg) possesses an essential singularity at the origin of the complex ctg- 



plane which agrees with the results of Refs. 15|-|17l|. 



Similarly, the sets of higher-order functions of the APT expansions, |2tfc(x)}, 
{v4fc(a;)}, . . . map onto the nonpower sets {w^(as)}, {wf (ctg)}, ... of the functions^: 



<(as) = ( L = - ) , w^{a,) = A, (- 
«;F(as) = Hfe(lnr2A2 = - 



(2.19) 



The adjacent elements satisfy simple differential relations 

<i(«s) = ^^^^, /? = E,M,D,..., (2.20) 



which follow from Eqs. fl2.16p at the one-loop level. Meanwhile, the same-order functions 
from different sets are interrelated with each other by the integral transformations following 
from Eq. (I2.10p . Besides, all the Euchdean functions («) possess an essential singularity 
at a = 0. 



^ Here Hfe(lnr^A^) denote the APT functions in the configuration representation which are related to Ak 



by the Fourier transformation 19 1 
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The regularity of the behavior of the timehke APT-invariant charge au and its effective 
"powers" 2tfc is provided by the summation [20^ of an infinite number of perturbative contri- 
butions related to the vr^-terms. In turn, the result of the 7r^-summation in the Minkowskian 
region, being extended into the Euclidean domain (12.101) . leads to the recovery of the non- 
analytic contributions of the form exp(— l/(/?o«/i)) which are "invisible" within the original 
perturbative expansion. Therefore, the Euclidean functions Ak contain both logarithmic 
and power terms in Q^. 

Apparently, the sets of the functions {w ^jas)}, appearing in the "ag-representation", 
are similar to Caprini-Fischer sequences (2l| which have been obtained proceeding from a 
different reasoning. 



E. Global APT and the "distorted" mirror 

In the studies of the hadronic processes at various energy scales one has to take into 
account the dependence of the theoretical results on the active quark flavor number Uf. 
Following the Bogoliubov method of massive renormgroup, the algorithm of the smooth 
matching was devised in Ref. This algorithm was employed for the analysis of the 



evolution of the running coupling in the range 3 GeV < Q < 100 GeV [23 . 

For the massless schemes, the matching of the invariant charge at the "Euclidean thresh- 
olds" is commonly employed . However, this procedure destroys the smoothness of the 
Euclidean function, and, therefore, violates its analyticity. 

The APT method opens a new opportunity for a self-consistent description of the observ- 
ables in domains corresponding to various numbers of active quarks nf |2^, j26|. Proceeding 
from the common matching condition, one defines the spectral functions p{a,nf) and em- 
ploys them in Eqs. (12. 6p and (12. lip afterwards. This procedure provides the analyticity of 
the Euclidean functions, whereas the Minkowskian functions turn out to be the piecewise 
smooth ones. Eventually, this results in the "global" APT-functions which "know" about 
all quark thresholds. 

Figure [3^, taken from Ref. [26|, depicts the global APT-charge in the Euclidean and 
Minkowskian domains. Figure [Hb demonstrates the "distorted mirror", i.e., that the Eu- 
clidean and Minkowskian functions possess a similar, but asymmetrical behavior. For com- 
parison with higher functions of the APT-expansions AkiQ"^) and 2tfe(s), the curves corre- 
sponding to the powers of the first APT functions are also given therein. 



F. Possible developments of the minimal APT 

As it has already been noted, the present version of the APT contains no additional 
parameters in comparison with the common RG-improved PT. We call it the minimal one. 
The straightforward application of the minimal APT in the low energy region (i.e., for the 
energies of the order of A) is not indisputable. It is worth noting that at such energies the 
effects due to the quark masses become considerable. In this situation, it is natural to modify 
the minimal APT by introducing new parameters. It is worthwhile to mention the so-called 
"synthetic" modification explored in Ref. j23]. Here the running coupling acquires the IR 
enhancement controlled by an additional parameter. In turn, this allows one to establish 
a link with the potential quark model. Other variants of the APT involve effective parton 
masses or modify the spectral representation by shifting the lower integration limit to the 
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FIG. 3: a) The global Euclidean and Minkowskian QCD running couplings for A^^^ = 350 MeV; b) 
the "distorted mirror" which elucidates the asymmetry of Euclidean and Minkowskian functions. 



two-pion threshold |l28j]. Besides, there are more formal ways of modification of the minimal 
APT at low energies [29|]. 

It is worthwhile to note that for numerical estimations one may employ the results of 
Ref. IJ] which also contains the expressions for the invariant expansion functions in terms 
of the Lambert function. Since these expansions are rather cumbersome, for practical appli- 
cations it proves to be convenient to employ simple approximate formulae. A simple explicit 
expression for the Euclidean QCD function was first proposed in Ref. 
one-parameter model 



30 in the form of a 



mod {Q ) — n 



1 

Jo 



+ 



exp(/2(a)) 



(2.21) 
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This model is based on the one-loop expression with the modified argument 
where a = 4. Recently this model has been essentially developed. It was shown sH, that 



h = l + ^Hl^ + a7r^), ^ = lni^, (2.22) 



the other choice of parameter a = 2 allows one to approximate both the Euclidean and 
Minkowskian three- loop APT functions Ak, 21^, = 1, 2, 3 within the accuracy sufficient for 
the description of all the experimental data above 1 GeV. The recurrent relations, modified 
in a proper way, lead to simple expressions of the "one-loop" form f l2.17p for higher functions. 
This enables one to employ the new model based on the substitution (12.221) in the analysis 
of contemporary experimental data without technical difficulties. This model has been used 



in the analysis of the inclusive T-decay in Ref. [31] • It was revealed therein that a weak 
point of the theoretical processing of rather precise experimental data is the choice of the 
scale. 



III. PHENOMENOLOGICAL APPLICATIONS 

The analytic approach has been successfully employed in studies of many hadronic pro- 
cesses. The literature devoted to the applications of the APT method is rather vast. Below, 
we mention some important results. 

A. Inclusive r-lepton decay 

The r lepton is the only lepton which is heavy enough to decay into hadrons. Experimen- 
tal data on the inclusive r lepton decay into hadrons possess good accuracy in comparison 
with those of other hadronic processes. These data constitute a "natural ground" for testing 
the low-energy QCD. 

The importance of the analyticity in the description of the r decay can be elucidated by 
the following example. The experimentally measurable quantity Rr is related to the life-time 
of the T lepton. The accuracy of its measurement is about 1 %. At the same time, Rr can 
be represented as the integral of the imaginary part of the correlation function 

«-^r^(i-^)\i + 2^)i-nnW. (3.1) 

The principal difficulty of the theoretical analysis of Rr is due to the fact that the integration 
range in Eq. (13. ip involves the low energy region. Meanwhile, the standard perturbation 
theory is not valid in the IR domain. Besides, if Imn(s) is parameterized by the power 
series in a^, the integral in Eq. (13. ip does not exist. This latter difficulty can be avoided by 
proceeding to the contour integral in the complex s-plane 



where D{—s) is the Adler function defined in Eq. (12. 8p . However, the original Eq. (13.10 can 
be represented in the form of Eq. (13. 2p only if the correlator satisfies the required analytic 
properties. This latter condition holds in the framework of the APT. The inclusive r lepton 



decay was studied in Refs. (32j, [33 . 
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B. e^e -annihilation into hadrons 

The pro cess of e"*'e~-annihilation into hadrons was examined in the framework of the APT 



in Ref. [3J]. The measurable quantity here is R{s) which is the ratio of the hadronic cross- 
section to the leptonic one. The theoretical analysis of R{s) entails certain complications. 
The APT results presented below correspond to the so-called "smeared" function -Ra(s) (ssj . 
The parameter A specifies a minimal "safe" distance from the cut in the complex s-plane 
which guarantees the absence of difficulties of the theoretical description of resonances. 
The function Ra{s) reads 

Ra{s) = ^ [U{s + lA) - U{s - tA)] , (3.3) 

with A being a finite parameter. By making use of the dispersion relation for 11 (g^), one 
can express the function (13. 3p in terms of the measurable ratio R{s): 



The function R{s) in the integrand can be approximated by the relevant experimental data 
at low and intermediate energies and by its perturbative prediction at high energies. This 
allows one to obtain the "experimental" curve for Ra{s). The integration of Eq. (13. 4p 
results in the smearing of the resonance structure of R{s). The other quantity, which we 
use in comparing the experimental and theoretical results, is the mentioned-above Euclidean 
Adler D-function. The "experimental" curve for -D(Q^) can be obtained in the same way as 
for Ra{s). 

In Figs, m and [5l the comparison of the APT results (denoted by the "AA" labels) with 
the experimental prediction is presented. Figure H] also shows the curve found with the PMS 
optimization (36l^of the third-order perturbative expansion. The experimental curve is also 
taken from Ref. |36|]. It is worth emphasizing that the APT method and PMS optimization 
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AA (3-loop) 
A A (1-loop) 
Experiment 



1 2 3 4 5 eg (GeV) 

FIG. 5: The "experimental" function Dj\{Q^) vs. the APT results. 



originate in different reasonings, though lead to close numerical results. Further we show 
that, contrary to perturbation theory, the APT results are practically independent of the 
subtraction scheme. Therefore, additional optimization of the scheme dependence in APT is 
not as important as in the perturbative case. The "experimental" curve for D/^ presented in 



Fig. [5] is taken from Ref. [37||. Figures [Hand [5] demonstrate good agreement of the theoretical 
APT results with the "smeared" phenomenological functions -Ra and Da which have been 
reconstructed by making use of experimental data on the e^e~-annihilation into hadrons. 



C. Renormalization scheme dependence 

An inevitable truncation of perturbative series leads to the well-known problem of the 
renormalization scheme dependence. It is worth noting that there is no firm criterion of the 
choice of the renormalization prescription. The partial sum of the perturbative series bears 
a dependence on the renormalization scheme, which is a source of the theoretical ambiguity 
in processing the data. In QCD, such ambiguity is the greater the smaller the energy scale 
is. Therefore, the analysis of the stability of the results should involve the investigation of 
both higher loop and scheme stability. 



The scheme stability of the APT results was first examined in Ref. 38j. In Fig. [6] (taken 
from Ref. [s^), the strong correction r(s) to the i?-ratio of e"^e~-annihilation into hadrons 
R{s) oc 1 + r(s) is shown. The curves presented therein were calculated within the APT and 
PT approaches at the three-loop level. In these calculations the widely-accepted MS-scheme 
and the so-called H-scheme were employed. The latter scheme is close to the former one in 
the sense of the cancelation index. The issue of the scheme dependence was discussed in 
detail in Ref. [soj. 

Figure M shows that for the energy of about few GeV the standard perturbative approach 
leads to large ambiguity due to the scheme dependence. The APT method drastically 
reduces the scheme dependence of the theoretical results. In particular, the APT curves 
corresponding to different schemes practically coincide. A similar result also takes place for 



the inclusive r lepton decay [33[ and for the sum rules of the deep inelastic lepton-hadron 
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FIG. 6: The strong correction to the i?-ratio within the analytic approach (AA) and the standard 
perturbation theory (PT) for various renormaUzation schemes. 



scatteriner 



D. Hadronic contribution to the muon anomalous magnetic moment and to the 
fine structure constant 



The APT method has recently been applied j40| to the description of the so-called R- 
related quantities. Among them the hadronic contribution to the muon anomalous magnetic 
moment plays an important role. The latter quantity (in the leading order in the electro- 
magnetic coupling a = aqEo) reads 



«r=o(z) / -K{s)Ris), (3.5) 

with K{s) being a known function. 

The hadronic contribution to the fine structure constant a can be represented in the form 



^«had(s) = -^^^/ —- — :• (3-6) 



The superscript "(5)" implies that the contributions of only first five quarks {u, d, s, c, 
and b) was retained here. 

Similar to the case of the r decay, the integration range in these expressions includes the 
low-energy region where the perturbation theory is inap plic able. Quantities (13. 5p and (13. 6p 



were evaluated in the framework of the APT in Ref. [40|. The latter also employs the 
assumption about the behavior of the quark mass function at low energy which is based on 
the nonperturbative solution of the Schwinger-Dyson equation. The obtained theoretical 
value 



al"^ = (698 ± 13) X 10"^° (3.7) 
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as{s) 




10 50 100 200 (GeV) 

FIG. 7: The analysis of the five-quark timehke domain. The difference between the points denoted 
by + and o, • is due to the TT'^-corrections. The solid APT curve corresponds to A^^^ = 270 MeV and 

as{Ml) = 0.124. The standard dot-dashed curve corresponds to = 213 MeV and as(M|) = 
0.118. 



is in good agreement with the phenomenological estimations of a^^'^ which employ the data 
on e^e~-annihilation and r decay [4l|, (i^ j. 

The hadronic contribution to the fine structure constant at the Z-boson scale evaluated 
in Ref. 

A«i5d(^l) = (278.2 ± 3.5) x 10"^ (3.8) 
agrees with the estimation of Ref. [i^l which uses the data on e"'"e~-annihilation into hadrons 



^4ad(^l) = (275.5 ± 1.9expt ± l.Srad) X 10"^ (3.9) 

Besides, the method of APT provides a reasonable description of some other quantities, 
e.g., the inclusive r decay characteristic in the vector channel i?^, and the functions D^IQ"^) 
and -Ra('S) corresponding to the r decay data. 



E. Some other applications 



In the framework of the APT the observables in the timelike domain (s-channel) can 
be represented as the nonpower expansion over the functions which retain the so-called tt^- 
terms. The analysis of the s-channel observables [1^ has revealed the following. For the 
energies above 50 GeV (rij = 5) the running coupling as gains the effective positive shift 
— +0.002 with respect to the standard two-loop (NLO) analysis. In the energy range 
10 -7- 50 GeV (n/ = 5) the value of this shift increases, namely, Aos — +0.003. This leads to 
a new value of the QCD invariant charge at the scale of the Z-boson mass: as(M|) = 0.124. 
The obtained results are presented in Fig. [7] taken form Ref. [26]. 

The convergence of the APT expansions is better than that of the perturbative series 
which is demonstrated in Table [TTl The results for the Gross-Llewellyn Smith sum rule for 
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the deep inelastic lepton-hadron scattering at \fQ^ ~ 1.76 GeV (Euclidean domain) and 
for the inclusive r lepton decay [Mr = 1.777 GeV, timelike region) are given therein. Other 
applications are presented in Table 2 of Ref. 261 . 



TABLE II: Various order contributions to observables within PT and APT methods 



Process 


Method 


1st order 


2nd order 


3rd order 


Gross-Llewellyn Smith 
sum rule (y^ ^1.76 GeV) 


PT 
APT 


65.1% 
75.7% 


24.4% 
20.7% 


10.5% 
3.6% 


Inclusive r lepton 
decay {Mr = 1.777 GeV) 


PT 
APT 


54.7% 
87.9% 


29.5% 
11.0% 


15.8% 
1.1% 



It is worthwhile to mention several other applications of APT. The analytic running 
coupling has been successfully employed in the analysis of the meson spectroscopy [43]. A 
recent (preliminary) result obtained by the Milano group shows that the form of the QCD 
interaction extracted from the light quarkonium spectrum as a function of below 1 GeV 
can be well approximated by the three-loop Euclidean , the value of the scale parameter 
A being close to its world average value^. 

The Euclidean running coupling and the APT nonpower expansion were employed in the 
description of the formfactor of the pion-photon transition with Sudakov suppression 44 
and the electromagnetic formfactors 45 



Besides, a strong reduction of the sensitivity of the 
results with respect to both the choice of the factorization scale 46| and the renormalization 
scheme 



47 was revealed. 



IV. FURTHER DEVELOPMENT OF THE APT: THE INELASTIC LEPTON- 
HADRON SCATTERING 

In the previous sections, the APT method has been applied in the analysis of the physical 
processes which can be described in terms of the two-point function, namely, the correlator 
of the quark currents. The analyticity condition was employed in the form of the Kallen- 
Lehmann representation. In this section, we address the inelastic lepton-hadron scattering 
which can be characterized by the structure functions of two scalar arguments. These 
functions possess rather complicated analytic properties. Nonetheless, the basic idea of the 
analytic approach turns out to be useful in this case, too. 



A. Jost Lehmann representation 

In the case of the inelastic lepton-hadron scattering, the general principles of the axiomatic 
QFT are accumulated in the Jost-Lehmann (JL) representation^ |48| for the structure func- 



tions. In the nucleon rest frame this representation reads [49 



W{u,Q^) = eiqo) j dudX^ S[ql - (Mu - - \^]ijiu, \^), (4.1) 



^ We thank Prof. G.Prosperi for providing us with this information. 

^ This representation is also known as the Jost-Lehmann-Dyson representation, see 
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the support of the distribution ^^(u, A^) being locahzed on the manifold 

P=|u|<l, ALn = M2(l-/r^)l (4.2) 

For the physical process of inelastic scattering the variables v and assume positive 
values. It is convenient to introduce a symmetric in v function, to be denoted by the 
same W{y^Q'^\ By making use of the radial symmetry of the distribution '?/'(u, A^) = 
iIj{p,X'^), which follows from the covariance, one can rewrite the JL representation in the 
covariant form [30| 

pl poo pi 

W{iy,Q^)= dpp^ dz5{Q^ + M^p" + A^ - 2zp^u^ + M^Q^ )^(p, A'). (4.3) 

Jo J\l,^ J-i 

In what follows we adopt the common notation = —q^ and v = qP, where P is the 
hadron momentum and q stands for the momentum transferred. 



B. Dispersion relation for the forward scattering amplitude 

Proceeding from the representation (14.31) . one can derive the //-dispersion relation (DR) 
for the virtual forward Compton scattering amplitude T(i/, Q"^) 

T{u,Q') = - t dpp" r dX' f dz ^{P^>^') ^ (4 4) 

In the complex z/^-plane this function has the cut along the positive semiaxis of real u"^. This 
cut starts at the point u^^^ which is determined by 

+ M-p- + A^ 



{\,p,z} 

that yields Umin = <5^/2. 

Thus, the DR at hand takes the form^ 



2zp 



1 f°° du'^ 

Q') = - I TY-T^—Wiuu Q')- (4.6) 



The DR for the i/-odd structure functions can be derived from the JL representation (14. ip 
in a similar way. 

In terms of the Bjorken variable x = Q'^/2i' DR ( 14. 6p takes the form 



^2^ _ rrr 2 dxi 1 



TT Jo Xi I - {Xi/X)^ 

Thus, T{x, Q"^) is an analytic function in the complex x-plane with the cut along the segment 
— 1 < X < 1. The DR (14.71) is suitable for establishing a link with the operator product 
expansion (OPE). 



This equation generalizes the DR for the real Compton effect; cf., e.g., with the results of Ref. [50 1. 
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A natural scaling variable for the JL representation reads SOj, 51 



1 cm^+upi 

In terms of this variable one can represent the dispersion integral (14. 7p in the following 
form:^ 

C. Analytic moments of structure functions 

For the JL representation a natural scaling variable has the form of Eq. (14. 8p . It can also 
be rewritten in terms of Bjorken variable x, namely, 



MniQ') = iQ'r-' I (4.12) 



, g^ + 4M^ ..... 

One can infer that for the physical processes the variable s assumes the values in the range 
between and 1. 

The variable s differs from both Bjorken and Nachtmann ones, which are commonly 
employed in the analysis of the deep inelastic scattering. However, it is this variable that 
leads to the moments which possess appropriate analytic properties. 

The s-moments of the structure functions can be defined as 

This function can also be rendered in the form of the Kallen-Lehmann representation 

which elucidates its analytic properties. The support of the weight function m„(o") is located 
on the semiaxis a > 0. The function m„(cr) can also be represented in terms of the original 
distribution ip{p, a) of the JL representation. 

The analytic moments (14. lip are related to the x-moments M^XQ"^) by j30| 

■^"^^^-r[(n + l)/2]^ kl M.+2.(Q). (4.13) 

In the asymptotic high energy region x-, s-, and ^-moments coincide with each other, since 
the power corrections l/((5^)" can be neglected. In the intermediate- and low-energy regions, 
where higher twist contributions become considerable, one has to distinguish between these 
moments. 

The moments of the structure functions (14. lip , being analytic functions in the complex 
Q^-plane with a cut, are convenient objects for employment within the analytic approach to 
QCD. The analyticity property is a consequence of the general principles of the local QFT. 







It is worth noting that in terms of another variable (e.g., the Nachtmann one ^ — 2x/(l 
a/1 + x'^AAP/Q'^)) such a representation does not exist. 
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D. Relation with the operator product expansion 

The identity of the structures of DR in x and s variables allows one to establish a relation 
of the analytic moments with the operator product expansion. The latter is commonly 
employed in the study of the Q^-evolution of the structure function moments. 

The x-moments of the structure function correspond to the case when only the Lorentz- 
structures of the form P,,, . . . P., are retained in the matrix element 

{P\d,,...,jP). (4.14) 

In this case, the OPE for the Compton amplitude leads to the expansion in powers of 
{qP)/Q^, i.e., in the inverse powers of x. The same expansion in the inverse powers of x can 
also be performed in the dispersion integral (14.71) . The relevant coefficients are determined 
by the x-moments. The comparison of these two expansions provides one with the relation 
of the x-moments with the OPE. 

In a general case, the symmetric matrix element (I4.14p contains the Lorentz-structures 
of the form 

{P^il ■ ■ ■ -Pju„}) M g^i^j{P^^ ■ ■ ■ -Pju„_2}; • • • • 

The ^-moments correspond to the choice of the operator basis, which involves the traceless 
tensors (i.e., such that the contraction of the metric g^^^j with (-P|0^i...^„ |-P) over any pair 
of indices vanishes) as the expansion elements. 

The dispersion relation (14. 9 p allows one to expand the Compton amplitude in inverse 
powers of s. If the OPE basis is chosen in such a way that an arbitrary contraction of the 
tensor (P|0^^...^„ |P) with the nucleon momentum P^. vanishes, then the OPE leads to a 
power series for the forward Compton scattering amplitude with the expansion parameter 
QfiQuiPfiPu — 9fj.uP'^){(l^)^^ which corresponds to expanding dispersion integral (14.90 in the 
inverse powers of s. This establishes a relation of the analytic moments of the structure 
functions with the OPE. It should be stressed that the orthogonality requirement of the 
symmetric tensor (P|0^^...^,JP) to the nucleon momentum P^^ unambiguously determines 
its Lorentz structure. 



E. Target mass effects 

The OPE method was applied to the description of the effects due to the mass of the target 



in Ref. [52|. This approach leads to the so-called ^-scaling, i.e., the parton distributions 
become functions of the ^-variable. The resulting expressions for the structure functions are 
inconsistent with the spectral condition. This situation reminds the ghost pole problem, 
namely, that an approximate solution contradicts the general principles of the theory. 

The solution of the spurious singularity problem was based on the Kallen-Lehmann rep- 
resentation. One can avoid the contradiction of the issue at hand with the spectral condition 
by accounting for the target mass corrections in the framework of the JL representation. 



This approach was implemented in Ref. j53| . Figure [H] depicts the ratio of the structure 



function P3, which retains the mass corrections, to the parton distribution at M'^/Q'^ = 1/2. 
The solid curve corresponds to the incorporation of the target mass effects by making use 
of the JL representation, whereas the dashed one corresponds to the ^-scaling method. One 
can infer from Fig. [8] that the structure function obtained within the ^-scaling method con- 
siderably deviates from that of the JL-method at large values of x. It is exactly the region 
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FIG. 8: Ratio of the structure function F3, which retains the massive corrections, to the parton 
distribution for MVQ^ = 1/2. The sohd curve corresponds to the incorporation of the effects due 
to the target mass by making use of the JL representation. The dashed curve corresponds to the 
^-scahng method. 



(a; ~ 1) where the ^-method comes into contradiction with the spectral condition. Since 
the accuracy of the experimental data is improving and such subtle effects as higher twists 
become important, in the theoretical studies one should originate in methods which are 
consistent with the general principles of QFT. 



V. CONCLUSION 

Renorminvariant PT is a basic tool for investigation of the QFT models as well as for 
practical calculations of the elementary particle interaction processes. The achievements 
of PT in QED and in the electroweak theory are well known. The perturbative compo- 
nent of QCD is substantial in the study of practically all the hadronic processes. In the 
low-energy region, where the nonperturbative effects become essential, the logarithmic con- 
tributions calculated within PT are usually supplemented with the power corrections of the 
nonperturbative origin. Thus, nonperturbative calculations involve PT as a component. As 
a result, PT affects the determination (from experimental data) of such nonperturbative 
quantities, as, for example, vacuum expectation values of quark and gluon condensates. 

The renormalized perturbative series contains large logarithms. Their presence keeps 
the terms of the series large at high energies, even when the expansion parameter (i.e., the 
coupling a^) is small. Nonetheless, since a perturbative expansion is not the final result of 
the theory, its certain modification is admissible. 

The RG method accumulates large logs within a new expansion parameter, namely, the 
invariant charge q;s(Q^). This results in perturbative expansions in powers of a^[Q'^\ The 
QCD invariant charge as(Q^) is small at high energies, that constitutes a crucial feature 
of QCD, namely, the asymptotic freedom (i.e., the quark-gluon interaction is suppressed 
at small distances). Thus, the principle of renorminvariance, which is implemented in the 
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RG method, allows one to modify the perturbativc series. The resulting expansions possess 
reasonable behavior in the UV domain (one has to keep in mind the asymptotic character 
of the perturbative series). These expansions are suitable for practical applications. 

However, at this stage the perturbative expansions entail the problems due to the unphys- 
ical singularities of the invariant charge. Such singularities contradict the general principles 
of the local QFT. The expansion in powers of asiQ"^) becomes unapplicablc at small en- 
ergies. Besides, at intermediate energies, which are important in QCD applications, the 
perturbative results contain large ambiguity due to the scheme dependence. 

The APT constitutes the next step in the modification of the perturbative expansion 
which eliminates the afore-mentioned difficulties. Here the principle of renorminvariance 
is supplemented with the causality principle. In the case of the quark current correlation 
function the latter principle is implemented in the analyticity condition which corresponds 
to the spectral Kallen-Lehmann representation. For the inelastic lepton-hadron scattering 
the general properties of structure functions are expressed in the integral Jost-Lehmann 
representation. 

In conclusion, we summarize the most important features of APT: 

• The effective Euchdean q;e(Q^) and Minkowskian aui-s) invariant charges are defined 
in a self-consistent way; they are free of unphysical singularities; they involve no additional 
parameters; the higher APT functions Ak and 21^ possess similar properties; 

• Euclidean invariant charge Q!e as a function of ag obeys an essential singularity 
exp(— 1//3oQ;s) at the origin; 

• It also has the infrared stable point q;e(0) = q;m(0) = l//9o which is independent of the 
scale parameter A; 

• In the framework of the APT the functional invariant expansions in powers of cis for 
observables are replaced by the nonpower expansions over the sets [Ak{Q^)} and {2lfc(s)}; 

• The dependence of the APT results on both the multi-loop (NNLO and N^LO) correc- 
tions and the choice of the renormalization scheme is drastically suppressed in comparison 
with that of the perturbative results. 

We conclude this paper with the utterance which was the epigraph of one of our papers: 

"Take care of Principles and the Principles will take care of you". 

We hope that this sentence is consonant with the creative work of Anatoly Alexeevich 
Logunov. 
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